Detection of BCS pairing in neutral Fermi fluids via Stokes scattering: 

the Hebel-SIichter effect 
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We consider the effects of superfluidity on the light scattering properties of a two component gas 
of fermionic atoms, demonstrating that the scattered intensities of the Stokes and anti-Stokes hues 
exhibit a large maximum below the critical temperature when the gas is superfluid. This effect, the 
light scattering analogue of the Hebel-Slichter effect in conventional superconductors, can be used 
to detect unambiguously the onset of superfluidity in an atomic gas in the BCS regime. 



Pacs Numbers: 03.75. Ss 05.30.Fk 32.80.-t 

Nearly a decade after the first experimental reports of 
atomic Bose-Einstein condensates (BEG) Jj, the trap- 
ping and cooling of gases with Fermi statistics has be- 
come one of the central areas of research within the field 
of ultracold atomic gases. Such gases offer the excit- 
ing prospect of examining the properties of interacting 
Fermi gases, including BGS superfluid states, with un- 
precedented accuracy and flexibility. 

Using Feshbach resonances, one can readily vary the in- 
teractions between atoms, allowing study of the crossover 
between a BEG of tightly bound molecules on the side 
of the resonance where < kpa <C 1, with a the scat- 
tering length and kp the Fermi momentum, to a BCS 
superfluid state on the other side of the resonance where 
< —kpo. <?C 1, with an interesting crossover regime 
in between |^, yl- Indeed, several groups have now re- 
ported clear experimental results for a BEG of diatomic 
molecules by looking at the momentum distribution of 
the gas Q . Experiments have also probed the a < side 
of the resonance using a fast magnetic sweep to the BEG 
side |3- Very recently, the a < side {—kpa ^ 3.3) was 
probed using a radio frequency (rf) transition to an un- 
paired hyperfine state |6|. The observed broadening of 
the rf spectral line was interpreted as arising from for- 
mation of Gooper pairs Q- For a trapped system this 
effect is complicated by the fact that the mean (Hartree) 
field also yields a significant broadening of the spectrum; 
the unpaired hyperfine state in general sees a different 
spatially dependent Hartree field than the paired ones, 
giving non-trivial energy shifts in the hyperfine transi- 
tions. The broadening of the rf line due to this effect 
has been observed experimentally Q, and it has been 
demonstrated for a spherical gas in the BGS regime that 
it introduces broadening of the rf spectrum at least com- 
parable to the effects coming from pairing Q- An inter- 
pretation of the broadening of the rf spectrum in terms 
of pairing effects only is therefore not straightforward for 
an inhomogeneous system. The absence of unambiguous 
signatures of the presence of superfluidity in the BGS 
limit (fci?|a| ^1), arises because the formation of large 
delocalized Gooper pairs does not affect the bulk proper- 
ties. Clear observation of BGS superfiuidity is one of the 
central problems for experimentalists. Suggestions to de- 



tect the onset of superfluidity include probing the collec- 
tive mode spectrum 10], the quantization of angula r mo- 
mentum [ij , off-resonance light scattering [13, [l3 j and 
by probing the dynamic structure factor using a scheme 
which avoids the complication due to the inhomogeneous 
Hartree field [3. 
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FIG. 1: Stokes (a) and anti-Stokes (b) scattering of an inci- 
dent laser beam. In (a) the laser, detuned off resonance by 
a frequency 5, excites a particle in the lower level | |) to an 
intermediate highly excited level | |), which de-excites, emit- 
ting a lower frequency photon, 7. The hyperfine splitting of 
the paired levels is ujm- In (b) the initial and final states are 
interchanged. 



One of the hallmark experimental tests of BGS the- 
ory in conventional superconductors was the observed en- 
hancement of the nuclear spin relaxation just below the 
transition temperature [l5j , an experiment which probed 
the detailed nature of the pair correlations. In this pa- 
per, we propose a related experiment to probe the pair 
correlations in trapped atomic gases, by looking at off- 
resonance inelastic (Stokes and anti-Stokes) light scatter- 
ing of a laser beam on a two component atomic Fermi gas. 
The present scheme differs from the earlier suggestions 
using off-resonant light scattering [13, UM crucially: we 
propose looking at Stokes and anti-Stokes inelastic light 
scattering involving a change of the hyperfine states of 
the atoms. The present proposal has several attractive 
features. First, there is no central coherent scattering 
peak, which is effectively insensitive to superfiuid corre- 
lations. Also, as we show, the effects due to superfluidity 
on the scattered spectrum can be strikingly large close 
to the critical temperature, Tc- Finally, since the present 



scheme simply flips the atoms between the two paired 
hyperfine states, it is not beset by compHcations due to 
non-trivial energy shifts coming from non-superfluid ef- 
fects, e.g., shifts in the Hartree energies 0jllJ|- 

We consider a trapped gas of atoms in two equally 
populated hyperfine states labelled by |cr), with a =t, i, 
interacting via an effective attractive interaction. We 
take the state | |) to have energy Whf (?l = 1 in this 
paper) above the state | t)j ^nd assume that below a 
transition temperature, T^,, the gas is BCS paired and 
superfluid. Consider a laser beam of frequency loi^ and 
wave vector k illuminating the gas. As illustrated in Fig. 
1, the light field can induce dipole radiation from the 
atoms by connecting the two hyperfine states a =t, i to 
an electronically excited state |e), which we take to have 
energy Wg above the state | t)- 

For large laser detuning, (5 = cj^ — Wg , the excited level 
|e) is not significantly populated and it can be adiabati- 
cally eliminated from the theory. The spectral intensity 
of the scattered light at position r is then 

(7 ■] (TO \J — OO >-' 

CT3CT4 

X « (riO)V'.3 (riO)V'i, (r2i)V'., (rs^)) (1) 

where the V'o- are the field operators for the two low-lying 
hyperfine states |cr) 'l7'|. Here /^^^|(r) is the scattered 
light intensity from a single atom, including the depen- 
dence of the atomic levels involved and the various di- 
rections of the experiment, and Ak is the change in mo- 
mentum of the scattered light compared to the incident 
light. The frequency lo — Log — cj^ is the difference be- 
tween the scattered light frequency los and the incident 
light frequency wl- 

Equation J^l describes two types of off-resonant light 
scattering processes. The first, elastic or coherent scat- 
tering is characterized by the initial and final atomic hy- 
perfine states being identical, corresponding to (T4 = (T3 
and Gi = o\. This process has been examined in detail 
by a number of authors and several effects of pairing have 
been identified. However, measurement of these effects 
is complicated by a large background coherent scattering 
intensit y w hich is largely independent of the state of the 
gas [Hill. 

The second type of scattering process is characterized 
by different initial and final atomic hyperfine states. If 
the initial state of the atom is | t) and the final is | J,), 
the emitted light frequency is reduced by tJhf from the 
incoming light frequency; this Stokes scattering process. 
Fig. la, corresponds to a\ = 0-4 =t and 02 = (T3 =i in 
Eq. (^. If, on the other hand, the initial state of the atom 
is I J,) and the final is | t) [ci = o'4 —[ and 02 = cr^ =t in 
Eq. QJ], the emitted light frequency is increased by Whf; 
this anti-Stokes scattering process is shown in Fig. lb. 
Stokes and anti-Stokes scattering therefore probe corre- 
lation functions that involve a hyperfine state "spin" flip 
[at positions ri and r2 in Eq.J^l]- As we show, the rate of 
these Stokes and anti-Stokes transitions involving atomic 



hyperfine state fiips are strongly affected close to Tc by 
BCS pairing of the states, in contrast to elastic light scat- 
tering which does not involve a hyperfine flip. We note 
that for such inelastic scattering to occur, dipole transi- 
tions between the electronically excited (orbital p) state 
|e) and the two (orbital s) hyperfine states should be al- 
lowed by the selection rules. That is, both the dipole 
matrix elements relevant for Stokes and anti-Stokes scat- 
tering, (e|d • £| t) and (e|d ■ £\ i), with d the atomic 
dipole operator, and £ the electric field of the incident 
laser, must be non-zero. 

From Eq. ^, the problem of calculating the 
scattered light intensity is reduced to evaluating 
the Fourier transform of the correlation function 
{tpl^ (riO)V'ff3 (riO)V42 {^2t)^Jai {r2t)) ■ Experiments on the 
pairing transition for ultracold Fermi atomic gases use a 
Feshbach resonance to enhance the atom-atom interac- 
tion, thereby increasing Tc- The gas is then best regarded 
as a molecular BEC on one side of the resonance and a 
weakly coupled BCS superfluid on the other side of the 
resonance, with a crossover region in between |2,I3- We 
are interested here in the problem of detecting the pres- 
ence of large delocalized Cooper pairs on the BCS side 
of the resonance, where it is adequate to use mean-field 
BCS theory. We expand the field operators in Bogoliubov 
eigenmodes [M,,(r),i;,, (r)] with energy i?,,, which can be 
obtained from a solution of the Bogoliubov-dc Gcnnes 
equations [ig. With this expansion, Eq. (^ for Stokes 
and anti-Stokes scattering yields, 

S{uj) ex Y^ [^\{u*^v*, - v*u*^,)au\^ 

nv' 
x(l-/„)(l-/„0<5(t^ + S„'+^„) 

+ \{ujju*^, + t;^u*/)AkP/,,(l - frj')SibJ + E^> - £;,,)|, 

(2) 

where /,, — [exp(/?£'^) + 1]^"'^ is the Fermi function and 
/ik — /d'^r exp(— ikr)ft,(r) denotes the Fourier transform. 
Equation Q yields the intensity of the Stokes line if cD = 
LO + ciJhf and that of the anti-Stokes line ii uj = uj — cjhf • 

In Eq. (O, we have, for notational simplicity, omitted 
the prefactor in S{u)) involving the atomic dipole matrix 
elements and the directional dependences of the exper- 
imental setup. The first two terms in Eq. 101 describe 
the creation and annihilation, respectively, of two quasi- 
particles. The third term describes the scattering and 
hyperfine "spin" fiip of a quasiparticle. For Stokes scat- 
tering, the quasiparticle is in the initial state {rj \) with 
energy E^ and it scatters into the state {rj' |) with en- 
ergy E' -fcjhf . Energy conservation for this process reads 
uol + E,f — ujs + E'^ + cjhf . For anti-Stokes scattering the 
initial and final quasiparticle states are (77 i) and (77' ]) 
with energies i?,, -f Whf and E' respectively. 

We first consider light scattering on a homogeneous 
system, which can be approximately realized experimen- 



tally for a trapped gas by focusing the laser beam on 
an area much smaller than the size of the atomic cloud. 
The quasiparticle eigenfunctions are then plane waves; 
for small frequency shifts, to <^ A, where A is the super- 
fluid gap, we can neglect the terms in Eq. (0) describing 
the creation and annihilation of two quasiparticles, and 
obtain 

S{uj) ex ^ |UqUq+Ak + Wgi'q+Akl^/gCl - /q+Ak) 

q 

xS{u; + UJM + Eq+Ak~ Eq), (3) 

with ul^{l + CjEq)/2 and v^ = {l-^ jEq)/2. Here 
^q — q^/'^rn — ^, and Eq = ,/^^_|_^2^ j^^ |.]^g limit. 



AA; <C kp, and Co <^ Akkp/m this expression reduces to. 
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where E' = E—uo. Equation Q is identical to the expres- 
sion for the nuclear relaxation rate for polarized nuclei 
in a conventional superconductor [l^- In particular, we 
note that the scattered intensity is logarithmically diver- 
gent for a) = in a infinite size superfluid system. This 
apparent divergence, which arises from the divergent den- 
sity of states at the Fermi surface for a superconductor, 
only appears in "spin" -flip processes such as Stokes and 
anti-Stokes scattering. In reality the divergence is medi- 
ated by finite size and finite lifetime effects. 

Had we instead considered elastic scattering processes 
where the initial and final hyperfine states of the atom are 
identical, the coherence factor in the last term in Eq. (0) 
describing quasiparticle scattering would read uu'* —vv'* . 
This factor vanishes at the Fermi surface for a) = 0, and 
the coherence term for elastic scattering exactly cancels 
the divergent density of states at the Fermi surface. For 
Stokes and anti-Stokes scattering on the other hand, the 
coherence factor in Eq. |(2Jl reads uu'* + vv'* ~ 1 and the 
divergence in the density of states at the Fermi surface 
due to pairing shows up in the response of the gas. The 
coherence factor refiects the fact that the coupling of the 
quasiparticles to a "spin fiip" perturbation is essentially 
the same in the superfiuid and normal phases. It is im- 
portant that the measurement be carried out a finite Ak 
such that Co ^ Akkp/m. For Ak = 0, the overlap in- 
tegrals in Eq. ^ simply yield a 6ri^ri' selection rule, and 
the scattered signal becomes oc ^ f^{l — f^), which does 
not exhibit a peak below T^ due to pairing. Measuring 
the intensity of the Stokes and anti-Stokes lines allows 
one to test the detailed nature of the pairing correlations 
reflected in the coherence factors, as was done for the 
BCS theory [l|. 

To illustrate the effect, we evaluate Eq. (0J numerically, 
with a pairing gap obtained by solving the BCS gap equa- 
tion as a function of T. We have chosen kplal = 0.3 and 
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FIG. 2: The scattered Stokes and anti-Stokes light intensity 
in a homogeneous system as a function of the reduced tem- 
perature T /T^, normalized to the scattered intensity at T^ . 
The inset shows the BCS gap A(r). 



a frequency shift Co/tp = 0.0001, where €p is the Fermi 
energy. Figure |21 shows the temperature dependence of 
the scattered light intensity of the Stokes and anti-Stokes 
lines. We see that close to Tc the intensity of light scat- 
tered from the superfluid state is signiflcantly larger than 
from the normal state. The intensity for scattering from 
the normal state is oc T . The scattered light from the 
superfluid is very large close to Tc due to the density of 
states effect described above, while it becomes exponen- 
tially suppressed for T — > 0, since the density of quasipar- 
ticles available for scattering scales as exp(— /3A). Due 
to the large peak in the scattering intensity below Tc, 
a Stokes-anti-Stokes experiment could clearly reveal the 
presence of pairing. Note that as Co decreases, the peak 
below Tc grows. 

To study the effect of the trapping potential, we now 
examine the scattered light from a gas in a spherical 
trap Vpot('') = muoTr^ l"^- The Cooper pairing is between 
atoms in time-reversed angular momentum states, (Ijrn) 
and (Z, —m), where / is the single particle orbital angular 
momentum, and m its component along the incident laser 
beam. We solve the Bogoliubov-de Gennes equations nu- 
merically using the method described in Ref. T9||. With 
the quasiparticle energies and wave functions obtained 
from this calculation, we then compute the scattered in- 
tensity from Eq. (j^J, where the quantum number 77 now 
stands for (n, I, m), with n the radial quantum number. 

In Fig. 2, we show the calculated intensities, for 
1.6 X 10"* particles trapped with a critical temperature 
kgTc ~ 0.09ep, u) = 0.003eF, and Ak = 2/lh with 
Ih — {ituot)^^^'^ the trap length. As in the homogeneous 
case, the scattered intensity from the superfluid gas has 
a large maximum below Tc- Again, this peak is due to 
the increased density of states close to the Fermi level in 
the superfluid phase for T close to Tc- 
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FIG. 3: The scattered Stokes and anti-Stokes light intensity 
as a function of T, in units of the scattered intensity at Tc, 
for a spherically trapped system. The inset shows the lowest 
quasiparticle energies as a function of angular momentum I. 



The shell structure of normal phase quasiparticle levels 
is less pronounced than in the non-interacting case since 
the Hartree field breaks the degeneracy of each harmonic 
oscillator level in the normal phase |23| . This introduces 
a dispersion in the quasiparticle energies as a function of 
angular momentum I. The pairing suppresses this effect 
bringing the quasiparticles closer in energy, and yielding 
a maximum in the density of states close to the Fermi 
level. This effect is shown in the inset in Fig. 3 where 
the lowest quasiparticle energies are plotted as a function 
of I for the superfluid state (x) and the normal state 
(+) for TjTc ^ 0.7. We see that the dispersion of the 



energy levels as a function of I is larger in the normal 
phase than in the superfluid phase. Again, it is crucial 
to choose a finite momentum shift Afc — 2/lh such that 
the total angular momentum I of a quasiparticle is not 
conserved in the light scattering process. However for 
m the angular momentum along the beam axis, a given 
quasiparticle with angular quantum numbers (Z, m) can 
scatter to any quasiparticle state with quantum numbers 
(r, m) with V 7^ I, and the peak in the density of states 
close to the Fermi energy in the superfluid state shows 
up in the intensity of the scattered light. For a very 
small momentum shift Afc ^ IJ^^ , the scattered signal 
would, as for the homogeneous case, from Eq. (0), be 
ex ^ fri(\ — fn), which does not exhibit any peak below 
Tc due to pairing. In contrast to the homogeneous case, 
the increase in the scattering intensity due to pairing is 
smooth at Tc due to the finite size of the spherical system. 
Again, the size of the peak increases with decreasing uj. 

In summary, we propose the detection of pairing in 
a two component atomic Fermi gas by looking at off- 
resonant light scattering. The onset of Cooper pairing 
yields a large peak in the intensity of the Stokes and 
anti-Stokes lines below Tc for both a homogeneous and a 
trapped gas. The proposed effect, which is the light scat- 
tering analog of the famous Hebel-Slichter effect on the 
nuclear relaxation rate in conventional superconductors, 
is unique in that the superfluidity gives an enhanced sig- 
nal close to Tc- The interpretation of the results it is not 
beset by non-trivial mean field energy shifts. 
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